On the dwarf galaxies rotation curves diversity problem 
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Antonino Del Popolo 


Abstract: In this paper, we show how baryonic physics can solve the problem of the striking diversity in dwarf 
galaxies rotation curves shapes. To this aim, we compare the distribution of galaxies of the SPARC sample, in 
the plane V2kpc-Viast (being V2kpc the galaxy rotation velocity at 2 kpc, and Viast that outermost one) with that of 
galaxies that we simulated taking account of baryonic effects. The scatter in the rotation curves in the V2kpe-Viast 
plane, and the trend of the SPARC sample’s, and our simulated galaxies’, distribution is in good agreement. The 
solution of the ’ diversity” problem lies in the ability of baryonic process to produce non self-similar haloes, contrary 
to DM-only simulations. We show also that baryonic effects can reproduce the rotation curves of galaxies like [C2574 
characterized by a slow rising with radius. A solution to the diversity problem can be obtained taking appropriately 
into account the baryon physics effects. 
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1 Introduction 


Among the predictions of the ACDM paradigm, agreeing very well with a plethora of observations [1—7]*, appears 
the forecast that dark matter (DM) haloes have a cuspy density profile, with pœ r~t, close to the halo center, given 
by the so called Navarro-Frenk-White (NFW) profile [10]. More recent findings confirm this result, however with a 
lower slope for the cusp [11]. The rotation curves (RCs) of dwarf, and low brighteness galaxies (LSB) are usually 
characterized by a more gentle increase than predicted by the NFW profile. This problem was noticed for the first 
time more than two decades ago [12-23]. 

From a more general point of view, the cusp/core problem is better defined in terms of the excess of DM predicted 
in the inner parts of the galaxies compared with the observed inner slope, and it can be connected to the Too-Big- 
To-Fail problem [24-26]. 

While dwarf galaxies usually have cored profiles, a more detailed study shows a significant spread in their RCs, 
and the existence of cuspy dwarfs. 

Despite disagreement on the above discrepancy, the inner profiles of several dwarf galaxies or LSBs are clearly 
not always flat [e.g. 27, 28], and [29] noticed in the THINGS sample a clear mass dependence of inner profiles. [30] 
first predicted such mass dependence. Inner slope dependence on the halo mass and on the ratio (Mjas+M.)/Mbotat; 
where M, is the stellar mass, was shown in more recent studies [31, 32]. 

In addition, currently, different techniques (e.g., spherical Jeans equation, multiple stellar populations techniques, 
Schwarzschild modeling) applied to those or similar objects sometimes give different results (e.g., the cuspy profile 
found by [33] in Fornax, while [34-36], found a core). Finally, discrepancies are evident for larger objects than those 
MW satellites (e.g., [37] found —0.17 < a < —0.01 in the case of NGC2976, while [38] found a = —0.90+0.15, and 
[39] found a = —0.53 +0.14 if stars were used to trace the potential, or a= —0.30 +0.18 if it was the gas). 
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From the studies discussed above, and several others, results show the existence of a range of profiles, and that 
no agreement on the exact dark matter slopes distribution can be reached based on morphologies, despite current 
improvements in kinematic maps. 

Recently, [28] quantified this diversity in dwarf galaxies RCs. They compared the circular velocity at 2 kpc Vorpe, 
given a fixed maximum in the circular velocity, Vmaz. The scatter in Vzkpc, in the Vmas range 50-250 km/s, spans a 
factor 3-4. 

Several authors proposed solutions to that problem, almost all relying on core formation process due to supernovae 
feedback (hereafter SNM). At this point we want to emphasize that despite the majority of studies converging on 
the idea that baryon physics leads to the formation of cores [17, 18, 31, 40-43], some studies arrive at the opposite 
conclusion [44-47]. While this disagreement can be due to different physical processes included in the simulations, it 
motivates one to be more careful in accepting simulation results. As for the diversity problem, it brings about one 
question: why wasn’t the diversity problem seen and solved by the large number of hydro-simulations run in the last 
decade, and especially in the past years, before [28] pointed it out, if the galaxies they formed were as realistic as 
claimed? Actually, some years before [28], the problem was discussed and solved by means of baryon physics, using 
a semi-analytical model in [17]. That core formation model, differently from the SNM, is related to the exchange 
of energy and angular momentum between gas clumps and DM through dynamical friction (dynamical friction from 
baryonic clumps model (DFBC))[17, 18, 31, 40, 48-57]. 

However, after the [28] analysis, several authors claimed the problem could be solved by the same simulations, 
based on the SNM, that were previously blind to it. In particular, [58] showed that similar simulations can solve the 
diversity problem with cores formation by baryons, while [28] could not. Even in the case that the core formation 
scenario through supernovae feedback can solve the diversity problem, such scenario encounters serious difficulties 
to explain the structure of objects like I1C2574, which displays a core extending to 8 kpc, where there is no star (the 
1CG2574 half-mass radius is 5 kpc). Problems like this prompted [59] to explore a solution based on self-interacting 
dark matter (SIDM): they found that the SIDM alone cannot solve the problem and thus somehow some baryonic 
physics must be introduced. 

In this paper, the distribution of galaxies produced by the GCDM in the Vaypc-Viast; Where Vast is the outermost 
radius, will be compared, following [28], to the SPARC data [60], a collection of high quality RCs of nearby galaxies. 

We expect the mass dependence of the inner structure of the galaxies, as shown in [31, 56], to give rise to a scatter 
in the Vogpe-Viast plane. Such scatter is not possible in the CDM scenario, producing self-similar DM haloes. 

In the first Sec. 2, the model and observations will be briefly detailed, and their confrontation 
discussed in Sec. 3, before concluding remarks in Sec. 4. 


2 Model and comparison with observations 


Our study of the diversity problem involved a subsample of the SPARC sample [60], which is a 
collection of nearby galaxies high-quality rotation curves, to determine V;,., and Vaste The subsample 
characteristics are described in Sec. 2.2. 

We then simulated 100 galaxies with the DFBC model [40, 57], with similar characteristics to 
our SPARC subsample, and with M, = 10” — 10" Mọ, in a ACDM cosmology according to the [61] 
parameters. 

Finally, we compared the SPARC subsample and simulated V2,,, Vast, as Summarized in the fol- 
lowing. 


2.1 Model 


The model simulating galaxy formation we used has been described in several papers [17, 18, 40, 
56, 57]. It provides a highly improved spherical collapse models from that described by [62-67], and 
includes the effects of dark energy [68-70], random angular momentum [e.g., 65, 67] produced by 
the random motions arising in the collapse phase, ordered angular momentum [e.g., 71—73] arising 
from tidal torques, adiabatic contraction [e.g., 74-77], gas and stellar clumps interactions with DM 
through dynamical friction [40, 48-55], gas cooling, star formation, photoionization, supernova, and 
AGN feedback [78—80, see the following]. 

It follows perturbations of diffuse gas (baryons) and DM, which will give rise to the proto-structure, 
from the linear to non-linear phase, through turn-around and collapse. Baryons fraction is set to the 
*universal baryon fraction” f, = 0.17 +0.01 [81] [0.167 in 2]. The baryons collapse compresses DM 


(adiabatic compression), steepening the DM profile [74, 75, 77]. 
If a DM particle is located at a given radius r <r; 


r,M,(r;) =r [M (r) + Mam(r)] “l 


[71, 74, 82] where M;(r;) is the initial dark halo distribution, then M, the final mass distribution of 
baryons (i.e. for example, an exponential disk for spirals, or a Herquist configuration ([83, 84] for 
elliptical galaxies), r the final radius, and Mam the final DM distribution, are obtained through solving 
iterativelly Eq. 1 [85]. This model can be improved to better reproduce numerical simulations by 
assuming conservation of the product of the radius by the inside mass for that orbit-averaged radius 
[75]. 

Radiative processes form baryons clumps, in turn collapsing towards the center of the galaxy and 
forming stars. Clumps formation, their life-time, and observation are discussed in [57]. 

Dynamical friction (DF) between baryons and DM transfers energy and angular momentum (AM) 
to the DM component [48, 49]. This gives rise to a predominantly outwards motion of DM particles, 
reducing the central density, and transforming the cusp into a core. [40, Appendix D] describes how 
DF is taken into account, inserting the DF force in the equation of motion [Eq. A14, 40] to affects 
structure formation. 

That mechanism, flattening the cusp, is amplified in the case of rotationally supported (spiral) 
galaxies through the ordered AM, L, acquired by the proto-structure through tidal interactions with 
neighbors [71, 86-88], and by random AM, j, originating by the random motions arising in the collapse 
phase [65]. 

The “ordered” AM is calculated from evaluating the torque 7(r), and integrating it over time ([71], 
equation 35; see also Sect. C2 of [40]). “Random” AM is taken into account by assigning a specific 
angular momentum at turnaround (for details see Appendix C2 of [40]). 

A classical cooling flow served as gas cooling mechanism |e.g., 89] [see Sect. 2.2.2 of 79]. The 
inclusion by [78] and [79, Sect. 2.2.2 and 2.2.3] of star formation, reionization and supernovae feedback 
were replicated. Following [79], reionization reduces the baryon fraction by 


fo 


Mir) = 2 
Piao (2 vir) [1+ 0.26M; (2) /Mo]? , ( ) 


with the virial mass M, and Mp is the ”filtering mass” [see 90]. We take the reionization redshift 
in the range 11.5-15. Our treatment of supernovae feedback also follows [91]. In that stage, each 
supernova explosion expels gas in successive events, lowering stellar density. The smallest clumps are 
destroyed by feedback soon after stars are formed from a small part of their mass [55]. 

AGN quenching becomes important for masses ~ 6 x 10''M, [92]. Its feedback was taken into 
account modifying [93] model as in Sec. 2.3 of [80], by forming a Super-Massive-Black-Hole (SMB) 
when the star density exceedes 2.4 x 10°M,/kpc*?, then accreting mass into it. 

Our model demonstrated its robustness in several ways: 


a. cusp flattening from DM heating by collapsing baryonic clumps is in agreement with previous 
studies [48, 49, 51—55]. 


b. galaxies density profiles correct shape [40, 94], this before the [41, 95] SPH simulations, and 
correct clusters density profiles [18] were predicted, and a series of correlations in clusters’ 
observations [96, 97] were reobtained [98]. 


c. inner slope dependence on halo mass [31], and on the total baryonic content to total mass ratio 
[18] were predicted, in agreement with [99]. In addition to these dependence, the inner slope 
was also found to depend on the angular momentum [18]. Finally, the correct DM profile inner 
slope dependence on the halo mass is explained over 6 order of magnitudes in halo mass, from 
dwarves to clusters [17, 18, 31, 40, 98], a range that no other model can achieve. 


2.2 Observational data 


The choice we made of galaxies observations was a subsample extracted from the Spitzer Pho- 
tometry & Accurate Rotation Curves (SPARC) [60]. Spanning large ranges in morphologies, surface 
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Fig. 1. Comparison of the prediction from our model with a selection of SPARC data, using (a) (top panel) the 
spherical circular velocity approximation; (b) (bottom panel) the disk plane cylindrical gravitational potential 
approximation. The full dots represent the SPARC sample, the open violet squares our simulated galaxies, the 


violet line represents the mean trend line, while the dashed line the expectation if a NFW profile described all 
haloes. 


brightnesses, and luminosities, and presenting new surface photometry at 3.6 um and high-quality 
rotation curves from previous HI/Ha studies, the entire sample contains 175 nearby disk galaxies. To 
minimizes the scatter in the baryonic Tully-Fisher relation, and following [100], we assumed a mass- 
to-light ratio Y. =0.5Mo/Lo. To avoid a strong dependence of V»,,. on the extrapolation method, we 
excluded from the full sample all galaxies which rotation curves didn’t provide pairs of data points 
bracketting 2 kpc. We also applied to the sample the conditions 


a. selecting stellar masses similar to our simulated galaxies, and 
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Fig. 2. Rotation curve of IC2574 compared with our results. The SPARC rotation curve of IC2574 is presented as 
dots with error bars. Our simulations result yields the grey solid line, including the RCs from its stars (dotted 
green line), and gas disk (dashed red line) components, themselves summed into the total baryonic mass (blue solid 
line). 


b. selecting galaxies with inclination > 45° which are the most reliable RC data. 


3 Results 


Once the observed sample and model simulated galaxies populations were determined, we compared 
their location in the plane Vzkpe — Viast+ 
In Fig. 1, we compared the SPARC data with the result of our model. The models interpretation of 


GM(<r) 
r j 


circular velocity are based on the 3D spherical circular velocity approximation (Vzirc-spherical = 
with r the 3D radius and M(< r) its enclosed mass) in Fig. la and calculated using the disk plane 
cylindrical gravitational potential (Veirc—potential = 4 RÈ y with R the cylindrical radius and © the disk 


potential, restricted in the galactic plane z = 0) in Fig. 1b. In both panels of Fig. 1, the full dots 
represent the SPARC sample. The open violet squares mark our simulated galaxies, the violet line 
represents the mean trend line, while the dashed line the expectation if a NFW profile described all 
haloes. The plot shows the region occupied by our galaxies distributions predicted by our model are 
in much better agreement with the SPARC galaxy distribution and its scatter than the output of [28]. 
The maximum spread is observed in the velocity range at 50 < Vass < 100km/s <, which corresponds, to 
M, values at which the action of feedback is maximum. 

Denoting as utliers PARC galaxies having V>,,,. outside the +30 range determined with respect 
to our model trend line, there are two outliers, namely: 1C2574, UCG05750. The circular-potential 
velocity definition employed in Fig. 1b lowers the trend line, especially in the maximum feedback 
region, and gives rise to a small increase in scatter. This reintegrates, on one hand, the outliers 
IC2574, and UGC05750. 

The successful reproduction of a similar distribution and scatter to SPARC galaxies by our model 
reflects the ability to model the total and stellar mass dependence of galaxies DM density profile 
and RCs, as shown in [56]: galaxies with M, ~ 10°M, are cored with a very flat profile and inner 
slope a~0, while galaxies with M, <10* tend to be more cuspy (larger a) as the feedback mechanism 
efficiency, producing the cores, is decreased for lack of baryon clumps. Conversely, galaxies with 
M, > 10° also tend to be cuspier, this time as deeper potential well from the galactic centers larger 
quantities of stars steepen the cusp and lowers the feedback mechanism efficiency. is reached for 
rotational velocities < 10km/s, It is important to notice that our galaxies follow the SPARC trend 
at Vas Z 150km/s, thanks to our accounting of AGN feedback that counteracts baryons cooling, and 
modifies star formation. In addition, tidal interaction in our model makes it more environment 
dependent than several hydro-dynamic simulations, in which galaxies are usually isolated. 

A recent paper of [101] shows similar results. We want to stress, that the method used in the 
two papers is totally different: simulations ([101]), vs a semi-analytic method (this paper). The 


explanation of the ’diversity problem” is related to interaction of baryon clumps with dark matter, 
and not due to Supernovae Feedback. Our results are dependent from environment (galaxies are not 
isolated), and the effect of AGN feedback, differently from [101]. Those effects change the distribution 
of galaxies in the V>,,.-Vias, plane. There are also differences in the SPARC sample used. In that of 
[101] some SPARC galaxies, outliers in their case, were not plotted by the authors. 

As previously discussed, IC2574, and UGC05750 are no more outliers, once calculating V2. in the 
galactic plane, as seen in the bottom panel of Fig. 1. In Fig. 2 we showed how our model predictions 
for the case of IC2574 RC. There, the SPARC RC of IC2574 (dots with error-bars) is compared with 
our simulation’s most similar galaxy RC (solid line). The plot also displays the contribution to the 
RC coming from stars, as the dotted green line, from the gas disk, the dashed red line, and from the 
total baryonic mass with the blue line. It shows a good agreement between our calculated RC and 
the SPARC’s. The case of UCG05750 was not plotted since the RC fit using the previous correction 
gives as good results as for 1C2574. 


4 Conclusions 


In this paper, we studied the problem of the diversity of RCs shapes in dwarf galaxies [28]. To this 
aim, we simulated 100 galaxies with similar characteristics to a subsample of the SPARC compilation, 
and compared the distributions of galaxies of that subsample in the Vzkxpe-Viast plane with the simulated 
galaxies’. The distributions scatter, and trend show good agreement between the SPARC compilation 
and our galaxies. However, two outliers are present. Determining the circular velocity in the galactic 
plane, instead of from the spherical symmetrical evaluation, we showed that the two galaxies are no 
longer outliers. We also showed how the prediction of our model is in agreement with the observations 


for one of the outliers, namely IC2574. 
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